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Abstract 

Numerical simulations are performed for the stick-slip motion in the Burridge- 
Knopoff model of one block. Calculated amplitude increases with the driving veloc- 
ity. We argue that this effect can be a criterion to distinguish between the stick-slip 
and pure slip oscillations. 
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1 Introduction 

Understanding and control of frictional properties of matter remains a chal- 
lenge for a wide spectrum of sciences [1] . In particular, physical effects which 
appear in the low velocity regime are of interest [2,3,4,5]. While details of 
rubbing surfaces depend on mechanical properties of the probing system [6], 
our understanding of these detail is to be built up from simple models. 

Here we concentrate on the simplest version of the Burridge-Knopoff model 
[7,8] of one block, connected with a spring to a driving mechanism moving 
with a constant velocity vq. The block moves on a substrate under forces of 
spring and friction. The core of the problem is the friction force F dependence 
of the velocity v. As it was recognized by multiple authors, the origin of the 
instability of the uniform motion is that F(v) decreases at least near v = 0. 
This remains true for a chain of any number of blocks [9] . 
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We note that a more detailed model of friction is known [10,2] includes the d 
ependence of the friction force on the age of contact. This approach, known 
as Rice-Ruina model from names of its authors, has some shortcomings; in 
its simplest form it does not reproduce stable orbits. An argument, somewhat 
lengthy, was presented by us in Ref. [11]. However, more strong theorem can 
be proved very briefly, basing on the Dulac criterion [12]. The proof is shown 
in the Appendix. Here we are interested in the stable orbits, then we work 
with the Burridge-Knopoff model. 

It is often assumed that the only alternative for the uniform motion is the 
stick-slip motion, where the block periodically switches from motion to rest 
and motion again. However, as it was recently demonstrated within an analyt- 
ical approximation, within some range of the driving velocity the stick phase 
disappears, but still the uniform motion is not stable [13]. The motion was 
termed 'pure-slip oscillations'. The stick-slip, pure slip and uniform motion 
are then three distinct possibilities, one succeeding another when the veloc- 
ity t>o increases. From the results of Ref. [13], a criterion can be derived to 
distinguish stick-slip and pure slip phases. The amplitude of the oscillations 
increases with the driving velocity vo for the stick-slip phase and it decreases 
for the pure slip phase. 

In a recent report [14], oscillations of motion of a PMMA polymer sample 
scratched by a conical diamond intender were measured as dependent on the 
driving velocity and normal load. The amplitude and period of oscillations 
are found to decrease with the driving velocity. A question appears, if the 
oscillations are of stick-slip or pure slip kind? The range of velocities, which 
are of order of /j,m/s, do not allow to observe the stick effect directly. We note 
that in another experiment with PMMA ([15], Fig. 5), the amplitude of the 
oscillations increases with the velocity. On the other hand, one could argue 
that the validity of analytical calculations in Ref. [13] is limited by the assumed 
shape of the function F(v), which was a polynomial function of third order. 

The aim of this paper is to provide new argument in this discussion by nu- 
merical calculations, for the shape of the function F(v) other than in Ref. 
[13]. After Ref. [9], we assume this function to be proportional to (1 + av 
Alternatively we add a linear term bv to the same function, to reproduce the 
transition from the stick-slip to the uniform motion. However, here we are 
interested only in the stick-slip phase. The goal is to check if the amplitude of 
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the oscillations increases or decreases with the driving velocity vq. 

In the next section we present the assumed equations of motion, the model 
parameters and the numerical results. Conclusions which can be derived from 
these results are presented in Section III. The text is closed by the above 
remarked Appendix on the Rice-Ruina model. 

2 The model and the results 

For two versions of the formula for the friction force Fj(i>), i — 1, 2, the equa- 
tion of motion is 

mx = —k(x — v t) — Fi(x) (1) 

where m and x are the block mass and the block position, k is the spring 
constant, v t is the position of the driving mechanism, Fi(v) = /i W/(l + av), 
F 2 — Fi + bv, and fi W = -Fj(O) is the static friction force. This equation is 
equivalent to a dimensionless equation 

u = vt — u pu (2) 

where r = y/±t, u = v = (3v , (3 = x = ^,c = a/(3 and p = b/(3. 
This version, controlled by two (F = Fi) or three (F = F 2 ) parameters, is used 
for the numerical calculations. The values of the parameters can be reproduced 
at least partially from the experimental data [14]. The normal load W = 196 
mN, the velocity v is from 5 to 1000 fim/s. From the fitting of the results 
to the time dependence of the spring force (Fig. 2 in Ref. [14]) we can get 
fj,o = 0.46, a = 160 s/m and k = 62 N/m. The sample mass is about 1 g, but 
it is certainly much less that the mass m of the movable part of the setup, 
which cannot be deduced from the content of Ref. [14]. It seems reasonable 
to evaluate its order of magnitude to be about 1000 g. In Fig. 1 we show the 
time dependence of the spring force. This is to be compared with Fig. 2 of Ref. 
[14]. The accordance is good. This plot proves only that the fit is successful. 

However, the calculated velocity dependence of the amplitude of the spring 
force is entirely different from the experimental data. Both the curve obtained 
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Fig. 1. Calculated time dependence of the spring force. The range of force is taken 
as the amplitude. The results are to be compared with Fig. 2 in Ref. [14]. 

for 6 = and the experimental points are shown in Fig. 2. Other calculated 
curves are the driving velocity dependence of the period of oscillations and 
the period and amplitude dependences on the normal load. These curves agree 
with the experimental data in the sense that they increase or decrease in the 
same way as in experiment. More detailed comparison is difficult, because we 
do not know the parameters: the mass m and the coefficients a and b. How- 
ever, the results of the simulations make clear that the calculated amplitude 
increases with the driving velocity. The only exception is when the mass is 
very small, less than 1 g, but even in this case the obtained dependence is 
very weak and incomparable with the experimental data. 

In the same Fig. 2 we show also the data for p = 2.38. This value is small 
enough to keep the range of the driving velocity in the stick-slip regime. How- 
ever, our qualitative conclusions remain unchanged. 

A comparison of our results with the experimental data of Ref. [15] can be only 
qualitative. The experimental setup is different from a simple block plus spring. 
However, some quantities are analogous to the previous case. For example, the 
spring constant is substituted by an elastic constant of the PMMA sample. 
On the other hand, numerical value of the coefficient /3 remains unknown. The 
value of a is about 250 s/m [15]. The coefficient b is set to zero. In Fig. 3 we 
show the amplitude dependence on the driving velocity for three values of f3: 
the same as in Ref. [14] (f3 = 42 s/m), 158 s/m and 418 s/m. The amplitude is 
in arbitrary units. For all investigated cases (three shown plus several others) 
the amplitude either remains constant or increases with vq. The increase of 
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the experimental results of Ref. [15] can be reproduced for (3 = 158 s/m. The 
vertical scale of each curve depends on the coefficient [IqW. 



3 Conclusions 



Summarizing the experimental data [14,15], we are faced with two experiments 
with opposite behaviour of the amplitude of the oscillations with the driving 
velocity. In the case of Ref. [14], the amplitude decreases about two times 
when the velocity increases by two orders of magnitude. In the case of Ref. 
[15], it increases about two times when the velocity increases by one order of 
magnitude. This qualitative difference of behaviour calls for an explanation 
with qualitatively different mechanisms. 

We postulate a solution that the dependence of the amplitude on the driving 
velocity, as calculated in Ref. [13], should be treated as a criterion to distin- 
guish between the stick-slip and pure slip oscillations. With this criterion, the 
data of Ref. [14] should be assigned to pure slip oscillations, and the data of 
Ref. [15] - to the stick-slip phase. 

Acknowledgements The authors are grateful to Tristan Baumberger, Piotr 
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Appendix 

The equations of motion in the Rice-Ruina model are 



^(vot-x) =fi + B\n^- + A\n^- (3) 



* = 1"^ ( 4 ) 
where x — v$t is the block position with respect to the driving mechanism, vq 
is the driving velocity, k is the spring constant, W is the normal force, /z is a 
reference value of the friction coefficient for steady sliding at some velocity V r . 
During sliding, the microcontacts are refreshed, on average, after a distance 
D . The state of these microcontacts is described by the variable $, which 
interpolates between the time of stick for the block sticked and D /v for 



5 



steady sliding. Finally, $ = A)/K and A, B are unitless material constants. 
We note that B > A in the experimental data [16]. 

The equations can be transformed to an autonomous form. Denoting (x — 
v t)/D = a, $/$ = 0, V r t/D = r, V/V r = u, exp(-// M) = 7, KD /W = 
k, we get dimensionless equations 

a = — uj + 70~^ exp[— — a] (5) 

B Hi 

= 1- a exp[_ q,] (6) 

where the time derivative is over r. To apply the Dulac criterion [12], we 
need an auxiliary function g(a,<f>); here the simplest g — 1 is sufficient. The 
divergence is 

exp[~a] (7) 

and it is different from zero except at the transition point from the uniform 
to the stick-slip motion. Then, there is no periodic orbits if 5 — A — k ^ 0. 

References 

[1] Socoliuc A, Bennewitz R, Gnecco E, Meyer E. Transition from stick-slip to 
continuous sliding in atomic friction: entering a new regime of ultralow friction. 
Phys Rev Lett 2004; 92 (13): 134301. 

[2] Heslot F, Baumberger T, Perrin B, Caroli B, Caroli C. Creep, stick-slip, and 
dry-friction dynamics: Experiments and a heuristic model. Phys Rev E 1994; 
49 (6): 4973-4988. 

[3] Braun O M, Dauxois T, Peyrard M. Friction in a thin commensurate contact. 
Phys Rev B 1997; 56 (8): 4987-4995. 

[4] Albert R, Pfeifer M A, Barabasi A-L, Schiffer P. Slow drag in a granular 
medium. Phys Rev Lett 1999; 82 (1): 205-208. 

[5] Krylov S Yu, Jinesh K B, Valk H, Dienwiebel M, Frenken J W M. Thermally 
induced suppression of friction at the atomic scale. Phys Rev E 2005; 71 (6): 
65101. 



6 



[6] Rozman M G, Urbakh M, Klafter J. Controlling chaotic frictional forces. Phys. 
Rev. E 1998; 57 (6): 7340-7343. 

[7] Burridge R, Knopoff L. Model and theoretical seismicity. Bull Sei Soc Am 1967; 
57 (3): 341-371. 

[8] Carlson J M, Langer J S. Properties of earthquakes generated by fault dynamics. 
Phys Rev Lett 1989; 62 (22): 2632-2635. 

[9] de Sousa Vieira M. Chaos in a simple spring-block system. Phys Lett A 1995; 
198 (5-6): 407-414. 

[10] Rice J R, Ruina A L. Stability of steady frictional slipping. J Appl Mech 1983; 
50: 343-349. 

[11] Szkutnik J, Kulakowski K. Transition from Uniform Motion to Stick-Slip in 
the Rice- Ruina Model of One and Two Blocks. 17-th Marian Smoluchowski 
Symposium on Statistical Physics, Fundamentals and Applications, Zakopane, 
Poland, September 4-9, 2004 ( |nlin.SI/0412068 ). 

[12] Glendinning P. Stability, Instability and Chaos: an Introduction to the Theory 
of Nonlinear Differential Equations. Cambridge; Cambridge UP 1994. 

[13] Thomsen J J, Fidlin A. Analytical approximations for stick-slip vibrations 
amplitude. J Non-Linear Mech 2003; 38 (3): 389-403. 

[14] Zhang S L, Valentine J M. Stick-slip and temperature effect on the scratching 
of materials. Tribology Letters 2002; 12 (4): 195-202. 

[15] Briscoe B J, Chateauminois A, Chiu J, Vickery V. Acoustic noise emission in a 
model PMMA/steel fretting contact. In : Dalmaz D, et al., editors. Tribology 
Research: From Model Experiments to Industrial Problem. Tribology Series, 
39. Elsevier, 2001, p. 673-681. 

[16] Baumberger T, Caroli C, Perrin B, Ronsin O. Nonlinear analysis of the stick- 
slip bifurcation in the creep-controlled regime of dry friction. Phys Rev E 1995; 
51 (5): 4005-4010. 



7 



100 



a. 
E 
< 




100 

v (nm/s) 



Fig. 2. Amplitude of the friction force as dependent on the driving velocity vq. 
Calculations are performed for b = (solid line) and b = 100 s/m (dotted line). 
Experimental points are taken from Fig. 4 in Ref. [14]. 
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Fig. 3. Amplitude of the friction force as dependent on the driving velocity vq. 
Calculations are performed for b = and three different values of (3, from 42 s/m 
(highest curve) till 418 s/m (lowest curve). Experimental points are taken from Fig. 
4 in Ref. [15]. 
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